HIGH ACCURACY SEMIDEFINITE PROGRAMMING BOUNDS 
FOR KISSING NUMBERS 



HANS D. MITTELMANN AND FRANK VALLENTIN 

Abstract. The kissing number in n-dimensional Euclidean space is the max- 
imal number of non-overlapping unit spheres which simultaneously can touch 
a central unit sphere. Bachoc and Vallentin developed a method to find upper 
bounds for the kissing number based on semidefinite programming. This paper 
is a report on high accuracy calculations of these upper bounds for n < 24. 
The bound for n = 16 implies a conjecture of Conway and Sloane: There is no 
16-dimensional periodic sphere packing with average theta series 

1 + 7680g^ + 4320g'' + 276480g'5 -I- 614409"^ H 



1. Introduction 

In geometry, the kissing number in n-dimensional Euclidean space is the maximal 
number of non-overlapping unit spheres which simultaneously can touch a central 
unit sphere. The kissing number is only known in dimensions n = 1,2,3,4,8,24, 
and there were many attempts to find good lower and upper bounds. We refer to 
Casselman [3] for the history of this problem and to Pfender, Ziegler [Tl], Elkies 
[7], and Conway, Sloane 6 for more background information on sphere packing 
problems. 

Bachoc and Vallentin [I] develop a method (Section [2] recalls it) to find upper 
bounds for the kissing number based on semidefinite programming. Table 1 in 
Section [3l the main contribution of this paper, gives the values, i.e. the first ten 
significant digits, of these upper bounds for all dimensions 3, . . . , 24. In all cases 
they are the best known upper bounds. Dimension 5 is the first dimension in which 
the kissing number is not known. With our computation we could limit the range of 
possible values from {40, . . . , 45} to {40, . . . , 44}. In Section|3]we show that the high 
accuracy computations for the upper bounds in dimension 4 result into a question 
about a possible approach to prove the uniqueness of the kissing configuration in 4 
dimensions. 

Although acquiring the data for the table is a purely computational task we think 
that providing this table is valuable for several reasons: The kissing number is an 
important constant in geometry and results can depend on good upper bounds for 
it. For instance in Section[5]we show that there is no periodic point set in dimension 
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16 with average theta series 

1 + 7680g^ + 4320(7^ + 276480^^ + 61440g^ + ■ • • 

This proves a conjecture of Conway and Sloane [6, Chapter 7, page 190]. Fur- 
thermore, the actual computation of the table was very challenging. Bachoc and 
Vallentin [T] gave results for dimension 3, . . . , 10. However, they report on numeri- 
cal difficulties which prevented them from extending their results. Now using new, 
more sophisticated high accuracy software and faster computers and more compu- 
tation time we could overcome some of the numerical difficulties. Section[3]contains 
details about the computations. 



2. Notation 

In this section we set up the notation which is needed for our computation. For 
more information we refer to 1 . For natural numbers d and n > 3 let Sd{n) be the 
optimal value of the minimization problem 



d 

mm " 

fc=i 

ai, . . . ,ad e M, ai, . . . ,ad > 0, 



a 

iin{l + ^afc + 6n + (Fo,^o"(l,l,l)) : 



bii,bi2, 622 G M, ( IJJJ IJ^^ ) is positive semidefinite, 

Fk e m(''+i-*^) ^ (d+i-k) ^ positive semidefinite, fc = 0, . . . , d, 

q,qi e R[7i], deg{p + pqi) < d, p,pi sums of squares, 

4 

r,ri, . . . ,r4 G K[u, v, t], deg(r -I- ^p^r^) < d, r, ri, . . . , r4 sums of squares, 

1=1 

d d 



1 + akPl\u) + 26i2 + 622 + 3 ^(ft , Sl{u, u, 1)} + q{u) + p{u)qi{u) = 0, 

k=l k=0 

d 4 

^'22 + Y{Fk,Sk{u,v,t)) +r{u,v,t) + Yp,{u,v,t)ri{u,v,t) = O}. 



fe=0 



Here is the normalized Jacobi polynomial of degree k with P^(l) = 1 and 
parameters ((n — 3)/2, (n — 3)/2). In general, Jacobi polynomials with parameters 
(a, /3) are orthogonal polynomials for the measure (l-u)"(H-M)'^dwon the interval 
[— 1, 1]. Before we can define the matrices we first define the entry (i,j) with 
hj^O of the (infinite) matrix Y^f' containing polynomials in the variables u,v,w 

by 



Then we get S'k by symmetrization: S'k = gSo-'^^fc"' where a runs through all 
permutations of the variables u,v,t which acts on the matrix coefficients in the 
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obvious way. The polynomials p, pi , . . . , p4 are given by 

p{u) = -{u+l){u + ll2), 

Pl{u,V,t) ^p{u), P2{u,V,t) ^ p{v), P3{u,V,t) =p{t), 

P4{u, v,t) — 1 + 2uvt — — ~ . 

By {A, B) we denote the inner product between symmetric matrices trace(y4i?). 

In [l] it is shown that this minimization problem is a semidefinite program and 
that every upper bound on Sd{n) provides an upper bound for the kissing number 
in dimension n. Clearly, the numbers Sdin) form a monotonic decreasing sequence 
in d. 



3. Bounds for kissing numbers 





best lower 


best upper bound 




SDP 


n 


bound known 


previously known 




bound 


3 


12 


12 


sii(3) 


= 12.42167009 . . . 






[16] Schiitte, v.d. Waerden, 1953 


Sl2(3) 
Sl3(3) 
Sl4(3) 


= 12.40203212 . . . 
= 12.39266509 . . . 
= 12.38180947 . . . 


4 


24 


24 


sii(4) 


= 24.10550859 . . . 






[11] Musin, 2008 


Sl2(4) 
Sl3(4) 
Sl4(4) 


= 24.09098111 . . . 
= 24.07519774 . . . 
= 24.06628391 . . . 


5 


40 


45 


sii(5) 


= 45.06107293 . . . 






[1] Bachoc, Vallentin, 2008 


Sl2(5) 
Sl3(5) 
S14(5) 


= 45.02353644 . . . 
= 45.00650838 . . . 
= 44-99899685 . . . 


6 


72 


78 


fill (6) 


= 78.58344077 . . . 






1^ Bachoc, Vallentin, 2008 


Sl2(6) 
Sl3(6) 
Sl4(6) 


= 78.35518719... 
= 78.29404232 . . . 
= 78.24061272 . . . 


7 


126 


135 


sii(7) 


= 134-8824614 . . . 






yj Bachoc, Vallentin, 2008 


S12(7) 
Sl3(7) 
Sl4(7) 


= 134.7319671 . . . 
= 134.5730609 . . . 
= 134.4488169 . . . 


8 


240 


240 

[12] Odlyzko, Sloane, 1979 
[9] Levenshtein, 1979 


sii(8) 


= 240.0000000 . . . 


9 


306 


366 


sii(9) 


= 365.3229274 . . . 






[1] Bachoc, Vallentin, 2008 


Sl2(9) 
S13(9) 
S14(9) 


= 364-7282746 . . . 
= 364.3980087 . . . 
= 364.0919287 . . . 


10 


500 


567 


sii(lO) 


= 558.1442813 . . . 






1^ Bachoc, Vallentin, 2008 


Sl2(10) 
Sl3(10) 
Sl4(10) 


= 556.2840736 . . . 
= 555.2399024 . . . 
= 554.5075418 . . . 


11 


582 


915 


Sll(ll) 


= 878.6158044 . . . 






[12] Odlyzko, Sloane, 1979 


Sl2(ll) 
S13(11) 
Sl4(ll) 


= 873.3790094 . . . 
= 871.9718533 . . . 
= 870-8831157... 
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12 


840 


1416 


Sll 






1364.683810 






[12] Odlyzko, Sloane, 1979 


Si2 


'12) 




1362.200297 . . . 








Sl3 


'12) 




1359.283834 








-514 


'12) 




1357.889300 . . . 


13 


1130 


2233 


,Si 1 






2089.116331 






[12] Odlyzko, Sloane, 1979 


Si2 


'13) 




2080.631518 . . . 








.Si 


'13) 




2073.074796 








Sl4 


'13) 




2069.587585 . . . 


14 


1582 


3492 


Sll 


'14) 




3224 950751 






[12] Odlyzko, Sloane, 1979 


-512 


'14) 




3202.448902 . . . 








513 


'14) 




3189.127644 . . . 








Sl4 


'14) 




3183.133169 


15 


2564 


5431 


-Sll 


'15) 




4949.650431 






[I2] Odlyzko, Sloane, 1979 


Sl2 


'15) 




4893.479446 . . . 








Si 


'15) 




4876.037229 . . . 








Sl4 


'15) 




4866.245659 . . . 


16 


4320 


8312 


Sll 


'16) 




7515.952644 






[13] Pfender, 2007 


Sl2 


'16) 




7432.720718 . . . 








Si3 


'16) 




7374.093742 








Sl4 


'16) 




7355.809036 . . . 


17 


5346 


12210 


Sll 


'17) 




11568 41674 






[13] Pfender, 2007 


Sl2 


'17) 




11333.84265 . . . 








Sl3 


'17) 




11128.26227 . . . 








Sl4 


'17) 




11072.37543 . . . 


18 


7398 


17877 


Sll 


'18) 




17473.48016 






[12] Odlyzko, Sloane, 1979 


Sl2 


'18) 




17034.32488 . . . 








Sl3 


'18) 




16686.28908 . . . 








Sl4 


'18) 




16572.26478 . . . 


19 


10668 


25900 


Si 1 


'19) 




26397 34794 






_3, Boyvalenkov, 1994 


Sl2 


'19) 




25636.98958 . . . 








Sl3 


'19) 




25029.87432 








Sl4 


'19) 




24812.30254. .. 


20 


17400 


37974 


Sll 


'20) 




39045.32761 






[12] Odlyzko, Sloane, 1979 


Sl2 


(20) 




37844.10380 . . . 








,Si 


'20) 




37067.18966 








Sl4 


'20) 




36764.40138 . . . 


21 


27720 


56851 


Sll 


'21) 




58087.03857 






[3] Boyvalenkov, 1994 


Sl2 


(21) 




56079.21685 . . . 








Sl3 


;2i) 


= 


55170.03449 . . . 








Sl4 


;2i) 




54584.76757 . . . 


22 


49896 


86537 


Sll 


;22) 




87209.06261 . . . 






[12] Odlyzko, Sloane, 1979 


Sl2 


;22) 




84922.09101 . . . 








Sl3 


;22) 




84117.92103 . . . 








Sl4 


;22) 




82340.08003 . . . 


23 


93150 


128095 


Sll 


;23) 




128360.7969 . . . 






[3] Boyvalenkov, 1994 


Sl2 


;23) 




127323.7095 . . . 








Sl3 


;23) 




125978.7655 . . . 








Sl4 


;23) 




124416.9796. .. 


24 


196560 


196560 
[la Odlyzko, Sloane, 1979 
[9] Levenshtein, 1979 


Sll 


;24) 




196560.0000 . . . 



Table 1. New upper bounds for the kissing number (best known values are underlined). 
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Finding the solution of the semidefinite program defined in Section [2] is a com- 
putational challenge. It turns out that the major obstacle is numerical instability 
and not the problem size. When d is fixed, then the size of the input matrices stays 
constant with n; when n is fixed, then it grows rather moderately with d. 

There is a number of available software packages for solving semidefinite pro- 
grams. Mittelmann compares many existing packages in [10] . For our purpose 
first order, gradient-based methods such as SDPLR are far too inaccurate, and sec- 
ond order, primal-dual interior point methods are more suitable. Here increasingly 
ill-conditioned linear systems have to be solved even if the underlying problem is 
well-conditioned. This happens in the final phase of the algorithm when one ap- 
proaches an optimal solution. Our problems are not well-conditioned and even 
the most robust solver SeDuMi which uses partial quadruple arithmetic in the final 
phase does not produce reliable results for d > 10. 

We thus had to fall back on the implementation SDPA-GMP [S] which is much 
slower but much more accurate than other software packages because it uses the 
GNU Multiple Precision Arithmetic Library. We worked with 200-300 binary digits 
and relative stopping criteria of 10~^°. The ten significant digits listed in the table 
are thus guaranteed to be correct. One problem was the convergence. Even with 
the control parameter settings recommended by the authors of SDPA-GMP for "slow 
but stable" computations, the algorithm failed to converge in several instances. 
However, we found parameter choices which worked for all cases: We varied the 
parameter lambdaStar between 100 and 10000 depending on the case while the 
other parameters could be chosen at or near the values recommended for "slow 
but stable" performance. 

The computations were done on Intel Core 2 platforms with one and two Quad 
processors. The computation time varied between five and ten weeks per case for 
d = 12. An accuracy of 128 bits in SDPA-GMP did yield sufficient accuracy but did 
not yield a reduction in computing time. 

After the computations for the cases d = 11 and d — 12 were finished new 128- 
bit versions (quadruple precision) of SDPA and CSDP became available; partly with 
our assistance. These new versions do not rely on the GNU Multiple Precision 
Arithmetic Library. So the computation time for the cases d = 13 and d — 14: were 
reasonable: from one week to two and a half weeks. 

4. Question about the optimality of the D4 root system 
Looking at the values S(i(4) in Table 1 one is led to the following question: 
Question 4.1. Is lim Srf(4) = 24? 

d — >oo 

If the answer to this question is yes (which at the moment appears unlikely 
because we computed Si5(4) = 23.06274835 . . . ), then it would have two noteworthy 
consequences about optimality properties of the root system D4. 

The root system D4 defines (up to orthogonal transformations) a point config- 
uration on the unit sphere = {x € R'^ : x ■ x = 1} consisting of 24 points; it is 
the same point configuration as the one coming from the vertices of the regular 24 
cell. This point configuration has the property that the spherical distance of every 
two distinct points is at least arccosl/2. Hence, these points can be the maximal 
24 touching points of unit spheres kissing the central unit sphere S^. 
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If linid^oo Sd(4) — 24, then this would prove that the root system D4 is the 
unique optimal point configuration of cardinality 24. Here optimality means that 
one cannot distribute 24 points on so that the minimal spherical distance be- 
tween two distinct points exceeds arccosl/2. Thus, the root system D4 would be 
characterized by its kissing property. This is generally believed to be true but so 
far no proof could be given. 

Another consequence would be that there is no universally optimal point con- 
figuration of 24 points in S*"^ as conjectured by Cohn, Conway, Elkies, Kumar [5]. 
Universally optimal point configurations minimize every absolutely monotonic po- 
tential function. The conjecture would follow if the answer to our question is yes: 
Every universally optimal point configuration is automatically optimal and Cohn, 
Conway, Elkies, Kumar 5J show that the root system D4 is not universally optimal. 

5. Nonexistence of a sphere packing 

Our new upper bound of 7355 for the kissing number in dimension 16 implies 
that there is no periodic point set in dimension 16 whose average theta series equals 

(1) 1 + 7680g3 -I- 4320^"^ + 276480(7^ + 61440^^ -I- • ■ • . 

This settles a conjecture of Conway and Sloane 6, Chapter 7, page 190]. In this 
section we explain this result. We refer to Conway, Sloane 6 , Elkies ffl, and to 
Bowert [2] for more information. 

An n-dimensional periodic point set A is a finite union of translates of an n- 
dimensional lattice, i.e. one can write A as A = (AZ" -I- wi) U . . . U (AZ" + vn), 
with vi, . . . e M", and A : K" M" is a linear isomorphism. The average theta 
series of A is 

N N 

0aW = ^EE E withg^e-. 

i—l j — 1 v^l/^ 

This is a holomorphic function defined on the complex upper half-plane. A holo- 
morphic function / which is defined on the complex upper half-plane, which is 
meromorphic for z — > zoo, and which satisfies the transformation laws 

/(-1/z) = z^!{z), and /(z + 2) = /(z) for all z £ C with > 0, 

is called a modular forms of weight 8 for the Hecke group G{2). The expression ^ 
defines the unique modular form of weight 8 for the Hecke group G{2) which starts 
off with 1 -|- Oq^ + Oq^ . It is also called an extremal modular form, see Scharlau and 
Schulze-Pillot [HI 

If there would be a 16-dimensional periodic point set whose average theta series 
coincides with ([l]) then this periodic point set would define the sphere centers of 
a sphere packing with extraordinary high density (see [BJ Chapter 7, page 190]). 
At the same time the existence of such a periodic point set would show that the 
kissing number in dimension 16 is at least 7680. This is not the case. 
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